High Precision determination of the n, K, D and D s decay constants from lattice QCD 
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We determine D and D s decay constants from lattice QCD with 2% errors, 4 times better than 
experiment and previous theory: f Ds = 241(3) MeV, f D = 207(4) MeV and fojfo = 1.164(11). 
We also obtain ftc/f^ = 1.189(7) and (/d < ,//d)/(/k' //it) = 0.979(11). Combining with experiment 
gives Vus=0. 2262(14) and V ca /V c d of 4.43(41). We use a highly improved quark discretisation on 
MILC gluon fields that include realistic sea quarks, fixing the u/d,s and c masses from the n, K, 
and rjc meson masses. This allows a stringent test against experiment for D and D s masses for the 
first time (to within 7 MeV). 



The annihilation to a W boson of the D s , Dd, 7r or K 
meson is a 'gold-plated' process with leptonic width (for 
meson P of quark content ah) given, up to a calculated 
electromagnetic correction factor [l|, Q, by: 
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V a b is from the Cabibbo-Kobayashi-Maskawa (CKM) ma- 
trix and fp, the decay constant, parameterizes the ampli- 
tude for W annihilation. If V a b is known from elsewhere 
an experimental value for T gives fp, to be compared 
to theory. If not, an accurate theoretical value for fp, 
combined with experiment, can yield a value for V a b- 

fp is defined from (0|a7 M 75&|P(p)) = fpp^ calcula- 
ble in lattice QCD to handle quark confinement, and 
with QED effects omitted. The experimental leptonic 
decay rates for K and 7r are known very accurately and 
D and D s less so, but with expected errors shortly of a 
few percent. Accurate predictions from lattice QCD can 
be made now, ahead of these results and comparison will 
then be a severe test of lattice QCD (and QCD itself). 
This has impact on the confidence we have in similar 
matrix elements being calculated in lattice QCD for B 
mesons that provide key unitarity triangle constraints. 

A major error in lattice QCD until recently was missing 
sea quarks from the gluon field configurations on which 
calculations were done, because of numerical expense. 
This has now been overcome. The MILC collaboration Q 
has made ensembles at several different values of the lat- 
tice spacing, a, that include sea u and d (taken to have 
the same mass) and s quarks with the u/d quark mass 
taking a range of values down to m s /10. The sea quarks 
are implemented in the improved staggered (asqtad) for- 
malism by use of the fourth root of the quark determi- 
nant. This procedure, although deemed 'ugly', appears 
to be a valid discretisation of QCD 0]. 

For fry a large error can arise from the inaccuracy of 
discretisations of QCD for c quarks. Discretisation errors 
are set by powers of the mass in lattice units, m c a, and 



Lattice/sea 
uoami, uoam £ 



16 d x 48 
0.0194, 0.0484 
0.0097, 0.0484 
20 3 x 64 
0.02, 0.05 
0.01, 0.05 
24 3 x 64 
0.005, 0.05 
28 3 x 96 
0.0124, 0.031 
0.0062, 0.031 



valence 

ami, am a , am c , 1 + e 



0.0264, 0.066, 0.85, 0.66 
0.0132, 0.066, 0.85, 0.66 

0.0278, 0.0525, 0.648, 0.79 
0.01365, 0.0546, 0.66, 0.79 

0.0067, 0.0537, 0.65, 0.79 

0.01635, 0.03635, 0.427, 0.885 
0.00705, 0.0366, 0.43, 0.885 



n/a 



2.129(11) 
2.133(11) 

2.650(8) 
2.610(12) 

2.632(13) 

3.711(13) 
3.684(12) 



TABLE I: MILC configurations and mass parameters used 
for this analysis. The 16 x 48 lattices are 'very coarse', the 
20 3 x 64 and 24 3 x 64, 'coarse' and the 28 3 x 96, 'fine'. The 
sea asqtad quark masses (I = u/d) are given in the MILC 
convention with uo the plaquette tadpole parameter. Note 
that the sea s quark masses on fine and coarse lattices are 
above the subsequently determined physical value We 
make a small correction (with 50% uncertainty) to our re- 
sults for fir,K,D,D 3 to allow for this, based on our studies of 
their sea quark mass dependence and MILC results in It 
has negligible effect on our final numbers and errors. The lat- 
tice spacing values in units of ri after 'smoothing' are in the 
rightmost column 0, [nj. The central column gives the HISQ 
valence u/d, s and c masses along with the coefficient of the 
Naik term, 1 + e, used for c quarks 0]. 



this is not negligible at typical values of a. However, 
m c a is not so large that it can easily be removed from 
the problem using nonrelativistic methods as is done for b 
quarks, for example The key then to obtaining small 
errors for c quarks is a highly improved rclativistic action 
on reasonably fine lattices (where m c a ps 1/2). 

The FNAL and MILC collaborations previously ob- 
tained a prediction for fp> of 201(17) MeV and for fp> s 
of 249(16) MeV Q using the 'clover' action for c quarks. 
The 6%-8% error comes largely from discretization er- 
rors in the clover action. Our action is improved to a 
higher order in the lattice spacing and this means that 
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our results are both more accurate for the D and D s and 
more accurate for charmonium, allowing additional pre- 
dictive power. In addition, our formalism has a partially 
conserved current so we do not have to renormalise the 
lattice /d to give a result for the continuous real world 
of experiment. We can then reduce the error on fp> to 
2% and the ratio of decay constants even further. 

We use the Highly Improved Staggered Quark (HISQ) 
action, developed |7|] from the asqtad action by reducing 
by a factor of 3 the 'taste-changing' discretisation errors. 
Other discretisation errors are also small since, in com- 
mon with asqtad, HISQ includes a 'Naik' term to cancel 
standard tree-level a 2 errors in the discretisation of the 
Dirac derivative. For c quarks the largest remaining dis- 
cretisation error comes from radiative and tree-level cor- 
rections to the Naik term and we remove these by tuning 
the coefficient of the Naik term to obtain a 'speed of light' 
of 1 in the meson dispersion relation. The hadron mass 
is then given accurately by its energy at zero momentum, 
unlike the case for the clover action. 

In 0] we tested the HISQ action extensively in the 
charmonium sector, fixing m c so that the mass of the 
'goldstone' r\ c meson agreed with experiment. We showed 
that remaining discretisation errors are very small, being 
suppressed by powers of the velocity of the c quark be- 
yond the formal expectation of a s {m c a) 2 and (m c a) 4 . 
The charm quark masses and Naik coefficients used for 
different 'very coarse', 'coarse' and 'fine' MILC ensembles 
are given in Table [IJ For the s and u/d valence quarks 
we also use the HISQ action with masses in Table [J] 

We calculate local two-point goldstone pseudoscalar 
correlators at zero momentum from a precessing random 
wall source Q and fit the average to: 
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T is the time length of the lattice and t runs from to 
T. i denotes 'ordinary' exponential terms and ip, 'os- 
cillating' terms from opposite parity states. Oscillating 
terms are significant for D q states because m c — m q is 
relatively large, but not for tt or K. We use a number of 
exponentials, i and ip, in the range 2-6 and loosely con- 
strain higher order exponentials by the use of Bayesian 
priors [llj. Constraining the D and D s radial excitation 
energies to be similar improves the errors on the D. In 
lattice units, Mp = M and fp is related to ao through 
the partially conserved axial current relation, which gives 
fp = (m a + mt,)\/2ao/MQ. The resulting fitting error for 
all states is less than 0.1% on Mq in lattice units and less 
than 0.5% on decay constants. Full details will be given 
in a longer paper \\2\ . 

To convert the results to physical units we use the scale 
determined by the MILC collaboration (Table [H in 
terms of the heavy quark potential parameter, r\. r\ja 
is determined with an error of less than 0.5% and allows 
results to be tracked accurately as a function of sea u/d 



quark mass and lattice spacing. At the end, however, 
there is a larger uncertainty from the physical value of 
t i . This is obtained from the T spectrum using the non- 
rclativistic QCD action for b quarks on the same MILC 



ensembles 
GeV. 



13|, giving n = 0.321(5) fm, 



0.615(10) 
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FIG. 1: Masses of the D+ and D 3 meson as a function of the 
u/d mass in units of the s mass at three values of the lattice 
spacing. The very coarse results are the top ones in each set, 
then coarse, then fine. The lines give the simultaneous chiral 
fits and the dashed line the continuum extrapolation as de- 
scribed in the text. Our final error bars, including the overall 
scale uncertainty, are given by the shaded bands. These are 
offset from the dashed lines by an estimate of electromagnetic, 
m u 7^ rrid and other systematic corrections to the masses. The 
experimental results are marked at the physical md/m a - 

Our results are obtained from u/d masses larger than 
approximately three times the average m u j d of the physi- 
cal u and d quark masses. We obtain physical answers by 
extrapolating our results to the correct m u /d using chiral 
perturbation theory. In addition we have systematic er- 
rors from the finite lattice spacing values used. Since our 
results are so accurate we can also fit them as a function 
of a to extrapolate to the physical a = limit. These 
two extrapolations are connected through the discreti- 
sation errors in the light quark action and one way to 
treat those is by modifying chiral perturbation theory 
to handle them explicitly [6j. A more general approach, 
that allows us to handle light and heavy quark discreti- 
sation errors together, is to perform a simultaneous fit 
for both chiral and continuum extrapolations allowing 
for expected functional forms in both with a Bayesian 
analysis ll| to constrain the coefficients. We tested this 
method by using it to analyze hundreds of different fake 
datasets, generated using formulas from staggered chi- 
ral perturbation theory [141 ] with random couplings. As 
expected, we found that roughly 70% of the time the 
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extrapolated results were within one standard deviation 
(computed using a Gaussian approximation to \ 2 ) of the 
exact result from the formula, verifying the validity of 
our approach and of our error estimates. 

We fit our results to the standard continuum chi- 
ral expansions through first order [l5| . augmented by 
second and third-order polynomial terms in x q = 
Bom g /8(7r/ w ) 2 , where Bq = m 2 /(m u + rud) to leading 
order in chiral perturbation theory. The polynomial cor- 
rections are required by the precision of our data [2(| . 
We include D* — D mass difference terms in the D/D s 
chiral expansion and take the DD*ir coupling to have 
the value inferred at leading order from the experimental 
D* width, allowing for a 30% error from higher order ef- 
fects. We correct for the finite volume of our lattice from 
chiral perturbation theory, although only has correc- 
tions larger than 0.5%. Our corrections agree within 30% 
with those in Ref. and we take a 50% uncertainty in 
the correction. We fit the couplings in the chiral expan- 
sions simultaneously to our ir and K masses and decay 
constants. We do the same for the masses and decay con- 
stants of the D and D s . Given the couplings, we tune 
m u /d and m s so that our formulas give the experimental 
values for and itik after correcting for the u/d mass 
difference and electromagnetic effects 0, [l7| ■ 

We find that finite a errors are 2-3.5 times smaller with 
the HISQ quark action than with the asqtad action, but 
still visible in our results. We combine the extrapolation 
to a = with the quark-mass extrapolation by adding 
a 2 dependence to our chiral formulas. We expect leading 
discretization errors of various types: a s a 2 and a 4 er- 
rors from conventional sources; and oi^a 2 , a^a 2 \og(x Uj d) 
and or 8 arx u ,d from residual taste-changing interactions 
among the valence and sea light quarks. We do not have 
sufficient data to distinguish between these different func- 
tional forms, but we include all of them (with appropriate 
priors for their coefficients) in our fits so that uncertain- 
ties in the functional dependence on a 2 are correctly re- 
flected in our final error analysis. The a 2 extrapolations 
are sufficiently small with HISQ (1% or less for ir and K 
from fine results to the continuum; 2% for D and D s ) 
that the associated uncertainties in our final results are 
typically less than 0.5%. The combined chiral and con- 
tinuum Bayesian fits have 45 parameters for D/D s and 
48 for ir/K with 28 data points for each fit [13]. 

Fig. Q] shows the masses of the D and D s as a function 
of u/d quark mass. To reduce uncertainties from the scale 
and from c quark mass tuning, the meson masses were 
obtained from m£> q — m ric /2 + m 71cexpt /2. The lines show 
our simultaneous chiral fits at each value of the lattice 
spacing and the dashed line the consequent extrapolation 
to a = 0. The shaded bands give our final results: rriD, 
= 1.962(6) GeV, m D = 1.868(7) GeV. Experimental re- 
sults are 1.968 GeV and 1.869 GeV respectively. We also 
obtain (2m£> s — m ric )/(2m£ > ~ rn Vc ) = 1.251(15), in excel- 
lent agreement with experiment, 1.260(2) [2j. This last 
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TABLE II: Error budget (in %) for our decay constants and 
mass ratio, where A x = 2mo x — m Vc . The errors are defined 
so that it is easy to see how improvement will reduce them, 
e.g. the statistical uncertainty is the outcome of our fit, so 
that quadrupling statistics will halve it. The a 2 and m u / d 
extrapolation errors are the pieces of the Bayesian error that 
depend upon the prior widths in those extrapolations. 'm s 
evolution' refers to the error in running the quark masses to 
the same scale from different a values for the chiral extrapola- 
tion. The ri uncertainty comes from the error in the physical 
value of ri and the finite volume uncertainty allows for a 50% 
error in our finite volume adjustments described in the text. 



quantity is a non-trivial test of lattice QCD, since we are 
accurately reproducing the difference in binding energies 
between a heavy-heavy state (the r\ c used to determine 
m c ) and a heavy-light state (the D and D s ). Table ITT1 
gives our complete error budget for this quantity. 

Fig. [5] similarly shows our results for decay constants 
on each ensemble with complete error budgets in Table HT1 
Ik and f n show very small discretisation effects and good 
agreement with experiment when V u d is taken from nu- 
clear (3 decay and V us from Ki 3 decays . We obtain f„ 
= 132(2) MeV and f K = 157(2) MeV. Alternatively our 
result for fx;/ fir (1-189(7)) can be used, with experimen- 
tal leptonic branching fractions @, llil. to g ive V us . Using 
the recent KLOE result for the K [l9|, l20( we obtain V us 
= 0.2262(13) (4) where the first error is theoretical and 
the second experimental. This agrees with, but improves 
on, the Ki 3 result. Then 1 - V 2 d - V 2 S - V 2 b = 0.0006(8), 
a precise test of CKM matrix first-row unitarity. 

fjj and fu s show larger discretisation effects but a 
more benign chiral extrapolation. Our final results are: 
f Ds = 241(3) MeV, f D = 207(4) MeV and f D Jf D = 
1.164(11). These results are 4-5 times more accurate 
than previous full lattice QCD results 0] and exist- 
ing experimental determinations. An interesting quan- 
tity is the double ratio {JdJ fo)/(fK / It is esti- 
mated to be close to 1 from low order chiral perturba- 
tion theory 2l|. We are able to make a strong quanti- 
tative statement with a value of 0.979(11). Equivalently 
{*d./*d)/<JkIU) = 1.005(10), where $ = fy/M. We 
also obtain (.fnJ.fn)/(.fnJfn) = 1.03(3) using our pre- 
vious result for the B ratio [22j. The B ratio dominates 
the error but improvement of this is underway. 

The results for /d and fo s obtained from the experi- 
mental leptonic branching rates coupled with CKM ma- 
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FIG. 2: D, D s , if and n decay constants on very coarse, 
coarse and fine ensembles, as a function of the u/d quark mass. 
The chiral fits are performed simultaneously with those of 
the corresponding meson masses, and the resulting continuum 
extrapolation curve is given by the dashed line. For n, K we 
have x 2 /dof = 0.2 and for D, D s , xV dof = 0.6, each for 28 
pieces of data. The shaded band gives our final result. At the 
left are experimental results from CLEO-c [23|, [25| (with the r 
decay result above the fj, decay result for D s ) and BaBar [24|] 
(D 3 only) and from the Particle Data Tables @] for K and 7T. 
For the K we have updated the result quoted by the PDG to 
be consistent with their quoted value of V us ■ 



trix elements determined from other processes (assuming 
V cs = V u d) are also given in Fig. They are: fo s = 
264(17) MeV for_u decay and 310(26) MeV for t decay 



ml 

25|. 



from CLEO-c [23j and 283(23) MeV from BaBar 
and f D = 223(17) MeV from CLEO-c for /z decay 
Using our results for fo s and /d 3 //d and the experi- 
mental values from CLEO-c (2j| for fi decay (for consis- 
tency between fu s and fn) we can directly determine: 
V cs = 1.07(1)(7) and V cs /V cd = 4.43(4) (41). The first 
error is theoretical and the second, and dominant one, 
experimental. The result for V cs improves on the direct 



determination of 0.96(9) given by the PDG 
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